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Abstract 



In this paper, the one-dimensional time-fractional diffusion- wave equa- 
tion with the fractional derivative of order a, 1 < a < 2 is revisited. 
This equation interpolates between the diffusion and the wave equa- 
tions that behave quite differently regarding their response to a local- 
ized disturbance: whereas the diffusion equation describes a process, 
• i-H where a disturbance spreads infinitely fast, the propagation speed of 

^ the disturbance is a constant for the wave equation. For the time- 

^ fractional diffusion-wave equation, the propagation speed of a distur- 

bance is infinite, but its fundamental solution possesses a maximum 
that disperses with a finite speed. In this paper, the fundamental 
solution of the Cauchy problem for the time-fractional diffusion- wave 
equation, its maximum location, maximum value, and other impor- 
tant characteristics are investigated in detail. To illustrate analytical 
formulas, results of numerical calculations and plots are presented. 
Numerical algorithms and programs used to produce plots are dis- 
cussed. 
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1 Introduction 



Evolution equations related to phenomena intermediate between diffusion 
and wave propagation have attracted attention of a number of researchers 
since 1980's. This kind of phenomena is known to occur in viscoelastic media 
that combine characteristics of solid-like materials that exhibit waves propa- 
gation and fluid-like materials that support diffusion processes. In particular, 
analysis and results presented in Pipkin (1986)] and Kreis and Pipkin (1986)] 
should be mentioned. Being unaware of an interpretation of evolution equa- 
tions by means of fractional calculus, these authors still could provide an 
interesting example of relevance of the intermediate phenomena for models 
in continuum mechanics. 

Nowadays it is well recognized that evolution equations can be interpreted 
as differential equations of fractional order at the time when some hereditary 
mechanisms of power-law type are present in diffusion or wave phenomena. 
This has been shown in Mainardi and Tomirotti (1997)| and more recently 
in |Mainardi (2010)] , where propagation of pulses in linear viscoelastic media 
governed by constitutive equations of fractional order has been revisited. 

For analysis of the evolution equations of the type mentioned above, meth- 
ods and tools of fractional calculus, integral transforms, and higher transcen- 



dental functions have been employed in the pioneering papers |Wyss (1986) 



Schneider and Wyss (1989)], |Fujita (1990)| , |Gorenflo and Rutman (1994) 
Kochubei (1989)[|Kochubei (1990)|, and in the book [Priiss (1993)|. We also 



ment ion the papers jMainardi (1994) Mainardi (1996a) Mainardi (1996b)] 



and Mainardi and Tomirotti (1995)] , where fundamental solutions of the 



evolution equations related to phenomena intermediate between diffusion and 
wave propagation have been expressed in terms of some auxiliary functions 
of the Wright type that sometimes are referred to as Mainardi functions, 
Podlubny (1999)] , |Gorenflo et al. (1999)] |Gorenflo et al. (2000)] . 



see i.e. 



These functions as well as some techniques and methods of fractional calcu- 
lus, integral transforms, and higher transcendental functions will be used in 
our analysis. 
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It is well known that diffusion and wave equations behave quite differently 
regarding their response to a localized disturbance: whereas the diffusion 
equation describes a process, where a disturbance spreads infinitely fast, 
the propagation speed of the disturbance is constant for the wave equation. 
In a certain sense, the time-fractional diffusion-wave equation interpolates 
between these two different responses. On the one hand, the support of the 
solution to this equation is not compact on the real line for each t > for a 
non-negative disturbance that is not identically equal to zero, i.e. its response 
to a localized disturbance spreads infinitely fast (see Fujita (1990)| ). On the 
other hand, the fundamental solution of the time-fractional diffusion-wave 
equation possesses a maximum that disperses with a finite speed similar to 
the behavior of the fundamental solution of the wave equation. The problem 
to describe location of maximum of the fundamental solution of the Cauchy 
problem for the one-dimensional time-fractional diffusion-wave equation of 
order a, 1 < a < 2 was considered for the first time in Fujita (1990) |. Fujita 
proved that the fundamental solution takes its maximum at the point = 
±c a t a / 2 for each t > 0, where c a > is a constant determined by a. Recently, 
another proof of this formula for the maximum location along with numerical 
results for the constant c a for 1 < a < 2 were presented in Povstenko (2008)] . 
In this paper, we provide an extension and consolidation of these results along 
with some new analytical formulas, numerical algorithms, and pictures. 

The rest of the paper is organized as follows: in the 2nd section, problem 
formulation and some analytical results are given. Here we revisit the re- 
sults of Fujita and Povstenko and give some new insights into the problem. 
Especially the role of the symmetry group of scaling transformations of the 
time-fractional diffusion-wave equation in the maximum propagation prob- 
lem is emphasized. We derive a new formula for the maximum value of the 
Green function for the Cauchy problem for the time-fractional diffusion-wave 
equation. A new characteristic of the time-fractional diffusion-wave equation 
- the product of the maximum location of its fundamental solution and its 
maximum value - is introduced. For a fixed value of a, 1 < a < 2, this 
product is a constant for all t > that depends only on a. The product is 
equal to zero for the diffusion equation and to infinity for the wave equation, 
whereas it is finite, positive, and laying between these extreme values for the 
time-fractional diffusion- wave equation that justifies the fact that the time- 
fractional diffusion-wave equation interpolates between the diffusion and the 
wave equations. The 3rd section is devoted to presentation of the numeri- 
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cal algorithms used to calculate the fundamental solution and its important 
characteristics including the location of its maximum, its propagation speed, 
and the maximum value. Results of numerical calculations and plots are 
presented and discussed in detail. 



2 Problem formulation and analytical results 

This section is devoted to the problem formulation and some important an- 
alytical results. In particular, several representations of the fundamental 
solution of the Cauchy problem for the time-fractional diffusion-wave equa- 
tion in form of series and integrals are given. These representations are used 
to derive explicit formulas for the maximum location, maximum value, and 
the propagation speed of the maximum point. Besides, we give a new proof 
of the fact that a response of the time-fractional diffusion-wave equation to 
a localized disturbance spreads infinitely fast like in the case of the diffusion 
equation. 



2.1 Problem formulation 

In this paper, we deal with the family of evolution equations obtained from 
the standard diffusion equation (or the D'Alembert wave equation) by re- 
placing the first-order (or the second-order) time derivative by a fractional 
derivative of order a with 1 < a < 2, namely 

d a u d 2 u 

W = dx*' ^ 

where x G IR, t G 1R + denote the space and time variables, respectively. 

In (OJ), u = u(x,t) represents the response field variable and the fractional 
derivative of order a, n — 1 < a < n, n G IN is defined in the Caputo sense: 

1 Vr)— (2) 



dt a T(n-a) Jo x ' dr r - 

where T denotes the Gamma function. For a = n, n G IN, the Caputo 
fractional derivatives is defined as the standard derivative of order n. 

In order to guarantee existence and uniqueness of a solution, we must add 
to some initial and boundary conditions. Denoting by f(x) , x G IR and 
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h(t) , t G H + sufficiently well-behaved functions, the Cauchy problem for 
the time-fractional diffusion-wave equation with 1 < a < 2 is formulated as 
follows: 

J u(x, 0) = f(x) , — oo < x < +oo ; , > 

U(^oo,t) = 0, t > 0. U 

If 1 < a < 2 , we must add to ^ the initial value of the first time derivative 
of the field variable, u t (x, 0) , since in this case the Caputo fractional deriva- 
tive is expressed in terms of the second order time derivative. To ensure 
continuous dependence of the solution with respect to the parameter a we 
agree to assume 

u t (x,0) = 0, for 1 < a < 2. 



In view of our subsequent analysis we find it convenient to set v := a/2, so 
that 1/2 < v < 1 for 1 < a < 2. 

For the Cauchy problem, we introduce the so-called Green function Q c (x, t; v\ 
which represents the respective fundamental solution, obtained when f(x) = 
5(x), 5 being the Dirac 5- function. As a consequence, the solution of the 
Cauchy problem is obtained by a space convolution according to 

/+oo 
G c {x - £,t;v) /(£) 
-oo 

It should be noted that Q c (x,t; u) = Q c (\x\,t; v) since the Green function of 
the Cauchy problem turns out to be an even function of x. This means that 
we can restrict our investigation of the function Q c to non-negative values 
x > 0. 

For the standard diffusion equation [y = 1/2) it is well known that 

G c (x,t; 1/2) :=^(M) = ^e-*V(4t). 

In the limiting case v = 1 we recover the standard wave equation, for which 
we get 

g c (x, t; 1) := g?(x, t)= l - [5{x -t) + S(x + t)] . 

In the case 1/2 < v < 1, the Green function Q c will be determined in the next 
subsection by using the technique of the Laplace and the Fourier transforms. 
The representations of the Green function Q c are of course not new and have 
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been discussed in Mainardi (1994)| - [Mainardi (2011)] to mention only a few 



of many papers devoted to this topic. 

In this paper, we are interested in investigation of some important charac- 
teristics of the Green function Q c including location of its maximum point, 
its propagation speed, and its maximum value. 

2.2 Representations of the Green function 

Following |Mainardi (1994)| - |Mainardi (2011)| , some representations of the 
Green function Q c in form of integrals and series are presented and discussed 
in this subsection. 

In Mainardi (1994)| , the Laplace and Fourier transforms technique was em- 



ployed to deduce the following representation for the Green function Q c for 

2uxQ c (x,t;u) = F v (r) = vrM u (r), (4) 



x > and | < v < 1: 



where 

r = x/t u > 

is the similarity variable and 

i r v 1 . a — ra v 

F v {r):=— / e a - ra da , M u (r) := — / da 

y ' 2m J Ha ' y ' 2ni J Ha a 1 -^ 

are the two auxiliary functions nowadays referred to in the literature of Frac- 
tional Calculus as the Mainardi functions, and Ha denotes the Hankel path 
properly defined for the representation of the inverse of the Gamma function. 

Let us note that the similarity variable r = xjf plays a very important 
role in our analysis of location of a maximum point of the Green function 
Q c . In its turn, the form of the similarity variable can be explained by the 
Lie group analysis of the time-fractional diffusion- wave equation Q. As 



shown in Buckwar and Luchko (1998)], Luchko and Gorenflo (1998)] and 



Gorenflo et al. (2000)| , symmetry groups of scaling transformations for the 



time- and space-fractional partial differential equations have been constructed. 
In particular, it has been proved in Buckwar and Luchko (1998)| that the 



only invariant of the symmetry group T\ of scaling transformations of the 
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time-fractional diffusion-wave equation (IT]) has the form r)(x , t, u) = x/t v that 
explains the form of the scaling variable. 

Using the well known representation of the Wright function, which reads (in 
our notation) for z G C 

Wx Jz):=— da=Y——i r , (5) 

where A > — 1 and fi > 0, we recognize that the auxiliary functions are 
related to the Wright function according to 

F„(z) = W- Vt0 {-z) = uzM v (z), M„(z) = W- Vtl -„(-z). (6) 

The formula ^ along with ^ provides us with the series representations of 
the Mainardi functions and thus of the Green function Q c (for x > 0): 

1 1 i oo t_ r jfv\n 

ft(«, %u) = — F.(r) = - M.(r) = - E ^JJ.^y (7) 

The formulas @-Q can be used to give a new proof of the known fact that 
the support of the Green function Q c is not compact on the real line for each 
t > 0, i.e. that a response of the time-fractional diffusion-wave equation 
with 1/2 < v < 1 to a localized disturbance spreads infinitely fast. Indeed, 
because the Wr ight function (|5|) is an analytical function for A > —1 and 
/i > (see e.g. |Gorenflo et al. (1999)| ) that is not identically equal to zero 
(Wx,fj,(0) = l/r(/i) > 0), the set of its zeros is discrete and has no finite limit 
points in the complex plane and thus on the real line. This means that the 
support of the function Q c {x,t; v) = W- Ut i- u (—x/t u ) is not compact on the 
real line for each t > 0. This fact was proved in |Fujita (1990)] using the 
representation ^ of Q c as a function depending on the similarity variable 
and the asymptotics of this function. 

Finally we mention another integral representation of the Green function Q c 
that can be found e.g. in Mainardi et al. (2001) | or Povstenko (2008)] : 



Q c (x,t;u) = — E 2u {—K 2 t 2u ) cos(xK,)dn, (8) 
where E a (z) is the Mittag-Lefner function defined by the series 

oo 
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The representation ([8]) can be easily obtained by transforming the Cauchy 
problem for the equation into the Laplace-Fourier domain using the 
known formula 

c\ d ^^ ] s\ = s a £{u(t) ] s}-]ru (k) {0 + )s a - 1 - k , n-Ka<n, (10) 



k=0 



with n 6 IN, for the Laplace transform of the Caputo fractional derivative. 
This formula together with the standard formulas for the Fourier transform 
of the second derivative and of the Dirac 5-function lead to the representation 



s 2u-l 



G c (k,s,u) = 2u ■ — , v = a/2 (11) 



+ K 2 



of the Laplace- Fourier transform Q c of the Green function Q c . Using the 
well-known Laplace transform formula (see e.g. Podlubny (1999)] ) 



L{E a (-t Q );s} 



g a-l 



and applying to the R.H.S of the formula (11) first the inverse Laplace trans- 
form and then the inverse Fourier transform we obtain the integral represen- 
tation ^ if we take into consideration the fact that the Green function of 
the Cauchy problem is an even function of x that follows from the formula 

m. 



2.3 Maximum points of the Green function Q c 

In Fig. 1, several plots of the Green function Q c {x; v) := Q c {x, 1; v) for differ- 
ent values of the parameter v [v = a/2) are presented (see the next section 
for description of numerical algorithms and programs used to calculate the 
numerical values of the Green function) . It can be seen that each Green func- 
tion has an only maximum and that location of the maximum point changes 
with the value of v. In Fig. 2, the Green function Q c (x, t; v) is plotted for 
v = 0.875 from different perspectives. The plots show that both the location 
of maximum and the maximum value depend on the time t > 0: whereas 
the maximum value decreases with the time (Fig. 2, right), the ^-coordinate 
of the maximum location becomes even larger (Fig. 2, left). Below we give 
a surprisingly easy explicit formula that connects the location of maximum 
and the maximum value for all times t > for a fixed u, 1/2 < v < 1. 
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Figure 1: Green function Q c (x;v) := Q c (x,l;u): Plots for several different 
values of v 
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Figure 2: Green function Q c (x,t; v): Plots for v = 0.875 from different per- 
spectives 



9 



The aim of this subsection is to present some analytical formulas that de- 
scribe both the location of the maximum of the Green function Q c , its value, 
and their interconnection, as well as the propagation speed of the maximum 
location. 

In the paper jFujita (1990)| , an elegant proof of the fact that the Green 
function Q c (x,t; v) of the Cauchy problem takes its maximum at the point 
x*{t,v) = ±c u t u for each t > 0, where c v > is a constant determined by 
u, 1/2 < v < 1, has been presented. His reasoning was as follows: Let us 
consider the Green function at the point t = 1 and for x > 0: Q c {x', v) '■= 
Q c {x,\\v). For 1/2 < v < 1 the function Q c {x\ v) is a stable probability 



density and the stable densities are all unimodal (see |Lukacs (1960)| ). This 
means that G c {%] v) takes its maximum at a certain point x** = c v with a 
constant c u depending on v. It follows from the formula Q that 

Qc(x; v) = l -M u {x) (12) 

and 

G c (x, t; v) = t -^M„(xt- v ) = t- v Q c {xt-»; v). (13) 

Because the function Gc(x; v) takes its maximum at the point the 
function Q c (x, t; v) has to take its maximum at the point that satisfies the 



relation xj,~ v = x^ = c v due to the formula (13). Thus the maximum point 
of the Green function Q c (x,t;i/) is moving with the time according to the 
formula 

z*(t) = cX, v = a/2. (14) 

As we see, the main argument in the Fujita's proof is dependence of the 
Green function from the similarity variable xt~ v '. 



This argument can be used to give an analytical proof of the relation (14) 



following an idea presented in Povstenko (2008)|. It should be noted that 



derivation of ( 14 ) given in Povstenko (2008) | is based on the integral repre 



sentation (|8j) and contains some divergent integrals that should be interpreted 
in one or another generalized sense. To avoid this, we present here another 
proof of (14) based on the representation ([7]) of the the Green function via 



the Mainardi function and not on the integral representation 

Because the Mainardi function M v is an analytical function for v < 1 as a 
particular case of the Wright function and because of (j7j), there exist partial 
derivatives of the Green function Q c (x, t; v) of arbitrary orders for t > 0, x > 
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and we can use the standard analytical method for finding its extremum 
points. We first fix a value t > and look for the critical points of the Green 
function Q c (x, t; v) that are determined as solutions to the equation 

^-Q c {x,t;v) = t -^-M' v {xT v ) = (15) 

or to the equation 

M' u {xr v ) = 0. (16) 
Now we try to find a function x* = x* (t) that gives a solution to the equation 



(16) for each t > 0. The equation (16) can be interpreted as an implicit 



function that determines the function {t). Let us find the time- 

derivative of X* clS db derivative of an implicit function: 



a l M' v {xT v ) -vt- v - x xM'l(xt- v ) 
Wx M> v {xt-) = t~ v M'lixt-v) " t 



dt J 



a 



x 
v- 



We thus obtained a simple differential equation for x*(t) with the solution 
x*(t) = Ct u , where C = = c v that is in accordance with the formula 



(14). Let us mention that the same method can be applied for any twice- 
differentiable function that depends on the similarity variable xt~ u or another 
one in form of product of the power functions in x and t. As it is known, the 
Green functions for many linear partial differential equations of fractional 
order possess this property and can be investigated by the method presented 
above. 

As mentioned in Fujita (1990)| , the maximum point of the Green function 



Gc(x, t; v) propagates for t > with a finite speed v(t, v) that is determined 
by 

v(t,v):=x',{t) = vc v t v ~ X . (17) 

This formula shows that for every u, 1/2 < v < 1 the propagation speed 
of the maximum point of the Green function Q c is a decreasing function in 
t that varies from +oo at time t = 0+ to zero as t — > +oo. For v = 1/2 
(diffusion) the propagation speed is equal to zero because of ci/2 = whereas 
for v — 1 (wave propagation) it remains constant and is equal to c\ — 1. In 
Fig. 3, some plots of the propagation speed of the maximum point of the 
Green function Q c are given for different values of v. For large values of t, the 
smaller the value of v is, the smaller is the propagation speed for the same 



time instant. Conversely, according to the formula (17), the smaller the value 
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1(T 4 1(T 2 10° 10 2 



t 

Figure 3: Propagation speed of the maximum point: Plot of v(t, v) for dif- 
ferent values of v in the log-lin scale 



of v is, the bigger is the propagation speed for the same time instant when 
t — > + . For example, the propagation speed for v = 0.505 becomes greater 
than the one for v = 0.55 for t < 3.04E — 24 (of course, this effect is not 
visible in the plot of Fig. 3). 

Now we determine the maximum value of Q c (x,t; v) in dependence of time. 
Let us denote the maximum value by Q*(t; v) and find it by using the integral 
representation (fsl) : 



G*(t;v) := G c (x*(t),t;u) = - E 2v (-kH 2v ) cos(c/k)4. (18) 

7T JO V ' 

The variables substitution r = t v K reduces the integral in (18) to the form 
g*(t; u) = — I™ E 2u (-T 2 ) cos(c,r) dr, (19) 



7T J 



i.e. the maximum value Q*(t; v) of the Green function can be written in the 
form 

g*(t;v) = m„r u , (20) 
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m 



V 



1 f 00 / \ 

- / E 2u (-T 2 ) cos(c u t) dr. (21) 

7T JO V ' 



It follows from the relations (14) and (20) that the product 



g*(t; v) ■ x*(t) =c v m v , < t < oo (22) 

is a constant that depends only on v or on the order a of the fractional 
derivative in the equation Q, i.e. that the maximum locations and the 
corresponding maximum values specify a certain hyperbola for a fixed value of 
a and for < t < oo. This fact is rather unexpected and one should look for 
its explanation from the physical and/or probabilistic viewpoints. Let us note 
that the product {?*(£; v) -x*(t) is equal to zero in the case v = 1/2 (diffusion 
equation) because the maximum point is always located at the point = 
and to infinity in the case v — 1 (wave equation) because the maximum value 
is always equal to infinity. The product values for 1/2 < v < 1 are finite 
and laying between these extreme values that justifies the fact that the time- 
fractional diffusion-wave equation interpolates between the diffusion and the 
wave equations. 

In Fig. 4, we give some plots of the parametric curve (x*(t), {?*(£; z/)) for 
< t < oo that is in fact a hyperbola for different values of v. The vertex 
of the hyperbola tends to the point (0,0) when v tends to 1/2 (diffusion 
equation) and to infinity when v — > 1 (wave equation). 



Another interesting and important curve is presented in Fig. 5, where the 
product c v m v of the maximum location and the maximum value of the Green 
function Q c [x]v) is plotted for 1/2 < v < 1. As we have seen above, the 
constants c v (maximum location of G c {x', vj) and m v (maximum value of 
Q c (x; v)) are decisive for the behavior of the Green function Q c (x, t]v) for 
all t > because the maximum locations and values of this function can be 
determined via these constants for any time point t > (see the formulas 
(14) and (20)). As we can see in Fig. 5, the product c v m v is a monotonically 
increasing function that takes values between (diffusion equation) and +oo 
(wave equation). For 0.56 < u < 0.99, the product varies between 0.1 and 
10 , i.e. it changes very slowly on this interval. For v — > 1/2 and v — > 1 the 
product c u m v goes to (diffusion equation) and to +oo (wave equation), 
respectively, very fast. 
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x.(t) 



Figure 4: Maximum locations and maximum values of Q c (x,t; v) for a fixed 
value of v. Plots of the parametric curve (x*(t), Q*(t;u)), < t < oo for 
different values of v in the lin-log scale 




0.5 0.6 0.7 0.8 0.9 

v 



Figure 5: Product of maximum locations and maximum values of Q c (x,t; v) 
for a fixed time t — 1: Plot of c v m u , 1/2 < v < 1 in the lin-log scale 
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3 Numerical algorithms and results 



In the previous section, some analytical results regarding the location of the 
maximum of the Green function Q c , its maximum value, and the propagation 
speed of the maximum point as well as the plots in Fig. 1 - Fig. 5 were 
presented. Because the analytical formulas derived in the previous section 
contain the constants c u and m u that we could not determine in analytical 
form, we used some numerical algorithms and MATLAB-programs for their 
calculation. These algorithms along with some numerical results and plots 
are presented in this section. 

To start with, we first discuss algorithms for numerical evaluation of the 
Green function Q c . Because Q c is a particular case of the Wright function 
(see formula ([6])), one can of course use the algorithms for the numerical 
evaluation of the Wright function suggested in Luchko (2008)] to evaluate 
the Green function Q c . 

Another approach to numerical evaluation of Q c we employed to produce our 
plots for this paper is to use the integral representation (J8j). To evaluate 
the Mittag-Leffler function E a in (j8f, we applied the algorithms suggested 
Gorenflo and Luchko(2002j] and the MATLAB ograms that implement 



m 



these algorithms and are available from Matlab File Exchange (2005)] . Be- 



cause the Mittag-Leffler function has for < a < 2 the asymptotics (see e.g 
Podlubny (1999)] ) 



E a {-x) = — r + 0{x~ 2 ), x -> +oo, (23) 

x 1 (1 — a) 

we can estimate the length of the finite integration interval in the improper 
integral ([!]) that allows to reach the desired accuracy e. Indeed, let A >> t~ v 
and 

1 2t~ 2v 
~eu\Y{l-2u)\ < ' 

Then the estimate 

\E,J-K 2 t 2u )\ < 



kH 2v |r(l - 2v) 



holds true for k > A because of the asymptotic expansion (23) and we have 



1 



7T 



I A 



Ei v ( —K 2 t 2v ) cos(xk) dn 



-2v 



2 1 2t 

K 'a 7tkH 2i '\T(1-2u)\ ~ A Trim -2u)\ < C ' 
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The integral 



TT 



J E 2v (—n 2 t 2u ^J cos(xk) dn 



with a finite value of A can be then evaluated using any of the known quadra- 
ture formulas. If A satisfies the conditions mentioned above, we get the 



estimate 



1 M / x 

— / E2 P [—K 2 t 2v ) cos(xk) cLk — G c (x, t; u) 

TT JO 



< e 



with the desired accuracy e that has been used for numerical evaluation of 
the Green function Q c . The results of the numerical evaluation of the Green 
function Q c for the time t = 1 and for different values of v are presented in 
Fig. 1. In Fig. 2, 3D-plots of the Green function are given for v = 0.875, 
1 <t<2, and < x < 3 that illustrate a typical behavior of Q c . 

As can be seen in Fig. 1 and as expected, {? c (x, 1; v) has a unique maximum 
for each 1/2 < v < 1 and the maximum location changes with v. Surpris- 
ingly, the maximum location does not always lay between zero (maximum 
location for the diffusion equation, v = 1/2) and one (maximum location 
for the wave equation, v = 1). Below we consider this phenomena in more 
details. 



As we have seen in the previous section (formula (14)), the location of the 
maximum point of the Green function Q c (x, t]v) depends on the constant 
c u , i.e. on the location of its maximum for t = 1. It is therefore very 
important to calculate c u numerically and to visualize the dependence of 
c v on u, 1 /2 < v < 1 . Because we already know how to calculate the Green 
function Q c (x, 1; v) and because it possesses a unique maximum point = c u , 
it is an easy task to find the maximum location e.g. with the MATLAB 
Optimization Toolbox. The results of the calculations are presented in Fig. 
6. 



Fig. 6 shows that the curve c u = c v {y) has a maximum located at the point 
v m 0.85. The value of the maximum is approximately equal to 1.28. It is 
interesting to note that for 0.69 < v < 1 the value of c v is greater than or 
equal to one. 

Finally we present results of numerical evaluation of the maximum value m v 
of the Green function Q c (x, t; u) at the time t — 1 as function of u, 1/2 < v < 1 
in Fig. 6. It follows from (20) that the constant m v determines the maximum 
value of Q c (x,t;u) at any time t > 0. For numerical calculation of m u , the 
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Figure 6: Maximum locations and maximum values of the Green function 
Q c (x; v): Plots of c v and m v for 1/2 < v < 1 



formula (21) was used. As expected, m u tends to infinity as v tends to 1 that 
corresponds to the case of the wave equation. 

Another interesting feature of the curve m v = m^fV) that can be seen in Fig. 
6 is that m v is first monotonically decreasing and then starts to increase. 
The minimum location of m v = m^iv) is at v m 0.61 and the minimum 
value is nearly equal to 0.25. Whereas m u changes very slow on the interval 
< v < 0.95, it starts to rapidly grow in a small neighborhood of the 
point v — 1. It should be noted that despite of the fact that the curves 
m u = m v {v) and c v = c u (u) are not monotone and possess a minimum and a 
maximum, respectively, the product c v m u is a monotone increasing function 
for all v, 1/2 < v < 1 (see Fig. 5). 

Finally we note that most of the analytical results and numerical algorithms 
presented in this paper can be extended and applied for more general par- 
tial differential equations of fractional order in two and three dimensions. 
Another important point that should be investigated in future would be to 
give a physical and/or probabilistic interpretation of the results regarding 
the maximum location and maximum value of the Green function Q c of the 
Cauchy problem for the time-fractional diffusion-wave equation. Especially, 
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an appropriate interpretation of the maximum of c v that is greater that 1 
(1 = Ci) and the minimum of m v that is smaller than TO1/2 ~ 0.28 would 
be desirable. The fact that the product of the maximum location x*(t) and 
the maximum value G*(t; v) of the Green function Q c (x,t; v) is equal to the 
constant c u m u for all t > should be properly interpreted, too. All these 
matters will be considered elsewhere. 
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